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ABSTRACT 

We present a study of the far-IR properties of a stellar mass selected sample of 1.5 < z < 3 
galaxies with log(M*/M ) > 9.5 drawn from the GOODS NICMOS Survey (GNS), the deep- 
est if -band Hubble Space Telescope survey of its type prior to the installation of WFC3. 
We use far-IR and sub-mm data from the PACS and SPIRE instruments on-board Herschel, 
taken from the PACS Evolutionary Probe (PEP) and Herschel Multi-Tiered Extragalactic Sur- 
vey (HerMES) key projects respectively. We find a total of 22 GNS galaxies, with median 
log(M*/M0) = 10.8 and z = 2.0, associated with 250 |Ltm sources detected with SNR > 3. 
We derive mean total IR luminosity logLiR(L0) = 12.36 ±0.05 and corresponding star for- 
mation rate SFRtr + uv = (280 ±40) Mq yr _1 for these objects, and find them to have mean 
dust temperature Jdust ~ 35 K. We find that the SFR derived from the far-IR photometry com- 
bined with UV-based estimates of unobscured SFR for these galaxies is on average more than 
a factor of 2 higher than the SFR derived from extinction corrected UV emission alone, al- 
though we note that the IR-based estimate is subject to substantial Malmquist bias. To mitigate 
the effect of this bias and extend our study to fainter fluxes, we perform a stacking analysis 
to measure the mean SFR in bins of stellar mass. We obtain detections at the 2 — 4a level at 
SPIRE wavelengths for samples with log(M*/MQ) > 10. In contrast to the Herschel detected 
GNS galaxies, we find that estimates of SFRtr + uv for the stacked samples are comparable 
to those derived from extinction corrected UV emission, although the uncertainties are large. 
We find evidence for an increasing fraction of dust obscured star formation with stellar mass, 
finding SFRtr/SFRtjv M^ 7±0 2 , which is likely a consequence of the mass-metallicity re- 
lation. 



2 M. Hilton et al 



1 INTRODUCTION 

Star formation rates (SFRs) in galaxies can be measured using 
many different methods (see e.g. |Kennicutt| 1998) . The most easily 
accessible tracer at high-redshift (z > 1) is rest-frame UV emission, 
which correlates with the number of young, massive stars and hence 
the global SFR of a galaxy. However, in dusty galaxies, this re- 
quires a significant correction due to absorption of UV photons by 
dust, which can be estimated using the correlation between the UV 
and far-infrared (IR) luminosity ratio (Ljr/Luv, where Ljr is con- 
ventionally defined over the wavelength range 8 — 1000 |^m) and 
the UV slope ((3; typically determined from a power law fit of the 
form fi oc X$ between 1500 and 2800 A), which has been measured 
from local starburst galaxies (e.g. |[Meurer et a l. 1999; Calzett fet al.| 
2000 ). Observations at far-IR wavelengths are generally thought to 
quantify the amount of obscured star formation more directly, as 
UV radiation associated with young stellar populations is absorbed 
by interstellar dust and re-emitted at far-IR wavelengths, and have 
revealed that much of the star formation activity that occurred at 
z > 1 is obscured (e.g. | Perez-Gonzalez et al.|2005| Le Floc'h et"aT] 
[2005l|Caputi et al.|2007||Magnelli et aL|2009||2011) . 

Observations over the last decade spanning a wide range in 
redshift and galaxy environments have shown that stellar mass is 
a key parameter for predicting the properties of a given galaxy. At 
low redshift (z < 0.1), the most massive galaxies tend to be red 
and located in denser environments than bluer, lower mass galax- 
ies (e.g. |Baldry et al.|2006| >. Although the colour-density relation 
weakens as redshift increases, a strong colour-mass relation is still 
seen at z ~ 2 (e.g. Grutzbauch et al. 2011). For galaxies which 
are actively forming stars, SFR is seen to be correlated with stellar 
mass up to z ~ 3 (e.g. Daddi et al. 2007 ; Magdis et al. 2010, Oliver 



et al.|2010[ [Bauer et al.|2011[|Karim et al.|2011[|Rodighiero'eTal'] 



201 1 ). Environment, while certainly important (as seen by the dom- 
inance of early type, passively evolving galaxies in clusters), seems 
to be more weakly correlated with other galaxy properties in com- 
parison to stellar mass, particularly at high redshift (e.g. Pen g et alT] 
|2010 , Grutzbau ch et al . |20 1 1 \ . This suggests that studies of the as- 
sembly of stellar mass, much of which occurs in obscured bursts of 
star formation, are crucial for developing our understanding of the 
galaxy formation process. 

In this paper we use far-IR photometry from the Herschel 
Space Observatory (Pilb ratt et aT]|2010| ) HerMES ( [Oliver et aT] 
|2012| ) and PEP ( [Lutz et al.||2011| > key projects to investigate ob- 
scured star formation in a stellar mass selected galaxy sample: the 
GOODS NICMOS Survey (GNS; |Conselice et al.|201l) . The GNS 
sample is selected in the //-band and is estimated to be complete 
for galaxies with stellar masses down to log(M*/MQ) = 9.5 at 
z < 3 ( [Grutzbauch et al.|201 1 [[Mortlock et al.|201 l||Conselice et al.| 
201 1 ). Bauer et al. ( 201 1 ) carried out a study of star formation ac- 
tivity in the GNS sample over the redshift range 1.5 < z < 3. This 
coincides with the peak of cosmic star formation activity as mea- 
sured in the UV (e.g. Bouwens et al. 2009); note however that in 
the IR a flat plateau in the SFR density is seen from 1 < z < 2 
(e.g.|Magnelli et al.|2011[|Bethermin et al.|2011| ). The |Bauer et~aT] 
{201 1) study primarily used rest-frame UV luminosity (corrected 
for extinction according to the UV slope) to estimate SFRs. In ad- 
dition, they estimated obscured SFRs for the « 20 per cent of their 
sample that were detected at 24 \xm using the MIPS instrument 
on board Spitzer, finding that the inferred total star formation rate 
(SFRir + uv) is on average 3.5 times larger than the SFR derived 
from the UV-slope extinction corrected UV flux (SFRuv,corr)- This 
factor of 3.5 may be overestimated, as several previous studies have 



shown that while 24 \xm flux densities can be reasonably extrapo- 
lated to measure Ljr (and hence SFRjr) for galaxies at z < 1.5, this 
is not the case at higher redshift (e.g. Papovich et al. 2007 ; Murphy 
|et al. 2009 2011), where Ljr as estimated from 24 |xm photome- 
try alone can be a factor ~ 5 higher than Ljr measured for the same 
sources when additional longer wavelength photometry is available 
to constrain the SED fits. The discrepancy is greater for Ultra Lu- 
minous Infrared Galaxies (ULIRGs, which have Lir > 10 12 Lq). 
Similar results have been reported in studies using Herschel data 
(e.g. |Elbaz et al.|2010|[20TT]|Nordon et al.|2010||20T2] ). 

Star formation in the massive (M* > 10 11 Mq) galaxies on 
which most of the GNS fields are centred (see Section \2.\\ has 
been investigated using far-IR data from the Balloon-borne Large- 
Aperture Submillimeter Telescope (BLAST; |Viero et al.|2012| ) and 
Herschel ( Cava et al. 2010 ), who found that disk-like galaxies (se- 
lected by the use of the Sersic index) have significantly higher SFRs 
than spheroidal-like galaxies. In this work we aim to improve the 
characterisation of obscured star formation as a function of stellar 
mass at 1 .5 < z < 3, using the combination of Herschel photometry 
and the wide stellar mass range spanned by the full GNS sample 
(log(M*/M )>9.5). 

The structure of this paper is as follows. In Section [2] we give 
a brief overview of the GNS and the Herschel data used in this 
work. We investigate the properties of the GNS galaxies detected 
at 250 pm using Herschel in Section]!] We extend the study to lower 
luminosity galaxies through a stacking analysis which is presented 
in Section]?] We present our conclusions in Section[5] 

We assume a cosmology with Q. m =0.3, £2 a = 0.7, and Ho = 
70 km s _1 Mpc -1 throughout. All values for star formation rates 
and stellar masses assume a Chabrier (2003 ) Initial Mass Function 
(IMF), unless noted otherwise. 



2 DATA 



2.1 Galaxy sample 



The galaxy sample used in this work is taken from the GNS fCon-| 
selice et al. 2011), which consists of 60 F160W (//-band) pointed 
observations in the GOODS fields ( Giavalisco et al. 2004 ) using the 
NICMOS instrument on-board the Hubble Space Telescope (HST). 
The footprint of the GNS overlaid on the SPIRE 250 |^m maps is 
shown in Fig. [I] Each GNS field is ~ 50 /r on a side, and covers 
the region around one or more massive galaxies (M* > 10 11 Mq) 
at 1.7 < z < 2.9, initially selected using a variety of colour selec- 
tion techniques: distant red galaxies (DRGs; Papovich et al. 2006 ), 
IRAC extremely red objects (IEROs; |Yan et al.||2004| ) and BzK 
galaxies (Daddi et al.|2007|). While this selection is not homoge- 
neous, Conselice et al. ( 201 1) shows that this combination of colour 
selection techniques leads to an almost complete sample of massive 
(M* > 10 11 Mq) galaxies: no single one of these colour selection 
methods selects more than 70 per cent of the massive galaxy popu- 
lation that would be selected in a photometric redshift survey, while 
a subsequent stellar mass selection in these fields based on photo- 
metric redshifts found an almost identical massive galaxy sample 
to the initial colour-based selection ( Cons elice et al.|20TT] ). 

In addition to providing high-resolution near-IR photometry 
of the massive galaxies targeted in each GNS pointing, the depth 
of the survey allows galaxies with much lower stellar masses to be 
detected: GNS is complete for galaxies with stellar masses down 
to log(M*/M ) = 9.5 at z < 3 ( [Grutzbauch etal.|2011 



Mortlock 



|et al.||20lT] ). The stellar mass measurements are described in de- 
tail in Conselice et al. (2011); briefly, a grid of Bruzual & Chariot 
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Figure 1. Footprint of the GOODS NICMOS Survey (red) overlaid on the HerMES 250 urn maps of the GOODS-North (left) and GOODS-South (right) 
fields. Each GNS pointing is in the direction of one or more M* > 10 11 Mq galaxies at 1.7 < z < 2.9, and is about 50" on a side. 



(2003) stellar population models, with exponentially declining star 
formation histories (x-models, with 0.01 < t(Gyr) < 10), spanning 
a wide range in metallicity (—2.25 < [Fe/H] < +0.56), were fitted 
to the BVizH photometry for each galaxy. 

In this paper we use a sample of 860 1 .5 < z < 3 galaxies with 
log(M*/M0) > 9.5 drawn from the GNS (the redshift range is cho- 
sen to match previous analyses of this catalogue presented in, e.g. 
|Bauer et al.|2"0TT||Grutzbauch et al.|20lT||Mortlock et al.|201 1| >. We 
include galaxies with both spectroscopic and photometric redshifts, 
using the former where possible. We do not cut galaxies with low 
photometric redshift probability (P, the % 2 probability outputted by 
HYPERZ, the code used to compute the GNS photometric redshifts; 
Bolzonella et al. 2000), because a comparison of the spectroscopic 
and photometric redshifts showed that the scatter of the residuals 
is similar regardless of the cut in P (a z = 0.04^Jvhen using only 
galaxies with P > 95 per cent, compared to a z = 0.06 using the full 
sample; see |Grutzbauch et al.|2011||Bauer et al.|20lT) . Note that 
450 galaxies in this sample have P > 95 per cent. 

To reduce contamination of the sample by AGN, we remove 
galaxies found within a 2" matching radius of X-ray sources listed 
in the 2 Msec Chandra catalogues of | Alexander et al.| (2003 
GOODS-North) and |Luo et al.] ( |2008l GOODS-South). These cat- 
alogues have flux limits of « 1.4 x 10 -16 ergs cm -2 s _1 in the 2- 
8 keV band, and are therefore deep enough to allow sources brighter 



than L 



'X(2-8keV) ' 



> 4 x 10 42 ergs s _1 to be detected at z ~ 2 (assum- 
ing a power law spectrum with r = 2). 

Later in this paper, we measure SFRs for GNS galaxies from 
the Herschel IR data and compare these with UV-based SFR mea- 
surements from |Bauer et al.| |2011) for the same galaxy sample. 
Here, we briefly summarise the method used to estimate these UV- 
based SFRs. 



Baue r et aL] ( |2011| > estimated unobscured UV SFRs from K- 
corrected ACS Z85o~band flux measurements, applying the SFRtjv- 
L2800 relation of Kennicutt (1998), where L28OO 1S the UV luminos- 
ity at 2800 A. These were corrected for obscuration by dust using 
the UV slope ((3) to estimate the amount of extinction, where (3 was 
measured from the 1600 A and 2800 A luminosities of the best fit- 
ting model SED for each galaxy. A similar methodology to |Calzetti| 
|et al.||2000| was used to convert (3 values into extinction estimates 
at 2800 A. The typical uncertainty on the UV-slope extinction cor- 
rected SFR estimates (SFRuv,corr) is ~ 30 per cent ( [Bauer et"aT 
|20TI) . 



2.2 Infrared data 

The Herschel photometry used in this work is taken from two 
key projects. The PACS ([Poglitsch et al. 2010) Evolutionary Probe 
(PEP; |Lutz et al.|20lT| ) provides 100 and 160 [im data covering both 
GOODS fields, as well as 70 \im coverage of GOODS-S. Simula- 
tions show that in GOODS-N, the flux limits at 80 per cent com- 
pleteness are 4.5 and 7.0 mJy at 100 and 160 \xm respectively, while 
in GOODS-S the corresponding limits are 1.5, 2.0, and 4.8 mJy at 
70, 100, and 160 [im. We also use 250, 350, and 500 yim SPIRE 
( Griffin et al. 2010) imaging data which was obtained as p art of 
the Herschel Multi-Tiered Extragalactic Survey (HerMES) 2 1 Oliver 



|et al.|2010||20T2] ). Unlike the PACS data, the SPIRE data are domi- 
nated by confusion noise from unresolved background sources. The 
calibration of the SPIRE instrument is described in Swinyard et al. 
< [20T0l >. 

Photometry was performed on all the Herschel maps, using 
prior positions derived from the MIPS 24 |im catalogue of Mag- 
|nelli et ah] ( |2009| > for source extraction. This 24 |^m catalogue is 
extracted from the GOODS-Legacy program observations (PI: M. 



1 G z is defined as the scatter in the photometric redshift residuals, i.e. 5z = 

(^spec — -Zphot ) / ( 1 H~ ^spec ) • 
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Figure 2. Postage stamp (10" x 10") NICMOS F160W (#-band) images of GNS galaxies detected in HerMES with SNR > 3 at 250 urn. The red cross in 
each postage stamp marks the position of the corresponding matched object in the HerMES/PEP catalogue, which is extracted using MIPS 24 \im priors. The 
green circle indicates the 2" matching radius used for cross-matching between the two catalogues. 



Dickinson), and reaches a 5a depth of about 30 juJy. Note that by 
requiring a 24 Lim detection for source extraction in the Herschel 
maps, a small fraction of sources will be missed at the GOODS 
depth (< 10 per cent; e.g. | Roseboom et al.||2010[ |Magdis et al.| 
|2011| |Bethermin et al.|2012| ). A blind extraction might be able to 
find such sources, at the expense of significantly noisier photom- 
etry due to source blending. Fluxes in the PACS maps were mea- 
sured by fitting scaled PSFs at each object position, as in Magnelli 
|et al.| ( [2009] >. In the case of the longer wavelength HerMES data, 
photometry was performed on all sources simultaneously, with the 
24 Lim catalogue being used to provide reliable deblending, using 
a slightly modified version of the method described in Roseboom 



et al. (201 0]). The changes to the method are described in |Rose-| 
boom et al.| ( |2012| >; briefly, a global (rather than local) background 



estimate was used in producing the catalogues used in this work, 
and a different (and faster) model selection algorithm was used in 
the fitting procedure. Using this deblending method, reliable fluxes 
can be extracted close to the formal ^4 — 5 mJy SPIRE confusion 
noise (measured after a 3a con f source cut, where o con f is the con- 
fusion noise measured without this cut, Nguyen et al. 2010). The 
24 Lim prior positional information reduces the impact of confusion 
noise, and so the approximate 3a limit for the SPIRE catalogue at 
250 Lim used in this work is ~ 9 mJy in both fields. We use this 



catalogue to investigate the properties of GNS galaxies detected at 
250 Lim in Section|3] 



In Section [4] we present a stacking analysis of GNS galaxies 
in bins of stellar mass, and we use data from other infrared sur- 
veys to broaden the wavelength coverage outside of the Herschel 
bands. In both the GOODS-N and GOODS-S fields we use Spitzer 
MIPS 24 Lim maps, taken from the Far Infrared Deep Extragalactic 
Legacy Survey (FIDEL DR2; PI: Mark Dickinson; for GOODS- 
S) and the GOODS-Spitzer survey (for GOODS-N). In addition, 
in GOODS-N we make use of the combined AzTEC/MAMBO 
1160 Lim map of |Penner et"aT] ( [20Tl] ), while in GOODS-S we use 
the 870 Lim LABOCA map from LESS jWeifi et al.|2009| >. To sim- 
plify the stacking analysis, the MIPS and PACS maps (in surface 
brightness units) are cross-correlated with the appropriate area nor- 
malised point spread function such that each pixel in the result- 
ing map represents the maximum likelihood flux density (in Jy) of 
an isolated point source at that position. For the publicly available 
AzTEC/MAMBO and LESS maps, this operation has already been 
performed. 
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3 PROPERTIES OF SPIRE DETECTED GNS GALAXIES 



3.1 Cross-matching 

We cross-match the GNS catalogue with the HerMES/PEP cata- 
logue using a simple 2" matching radius. Since the HerMES/PEP 
catalogue was extracted using MIPS 24 \im prior positions, a small 
matching radius, appropriate to the astrometric accuracy achiev- 
able with MIPS at 24 \im, can be used (e.g. |Bai et al.|2007) . We 
select robust detections at 250 \im from the catalogue using cuts 
of ^250 > 3 x AS250> where AS250 is the flux uncer taint y (including 
confusion noise), i.e. S250 > 8 — 9 mJy (see Section]^]), and % 2 < 5 
(i.e. the goodness of fit of the source solution within the neighbour- 
hood of the source, see Roseb oom et al.|20"T0] >. We find that a to- 
tal of 22 GNS galaxies with 1.5 < z < 3 and log(Af*/M ) > 9.5 
are matched across both the GOODS-N and GOODS-S fields; this 
corresponds to « 2.5 per cent of the GNS sample within these stel- 
lar mass and redshift cuts. We note that if we repeat the selection 
at 350 |^m, we obtain a sample of 14 objects, only 1 of which is 
not in common with the 250 \im selected sample. This additional 
source is ID 283 in the GNS catalog, and has photometric redshift 
z p = 1.55±0.15 and stellar mass log(M*/M ) « 10.6. 

Fig.[5]shows 10" x 10" NICMOS F160W postage stamp im- 
ages centred on each detected GNS galaxy, with the position of the 
HerMES source and the 2" matching radius indicated. In almost all 
cases each GNS galaxy is unambiguously identified with the Her- 
MES source; there are only two cases (IDs 4180 and 5310) where 
two galaxies of similar brightness are located within the matching 
circle. We estimated the fraction of potentially spurious matches 
by randomising the positions of the sub-mm sources and repeating 
the cross matching procedure 1000 times. We found a mean num- 
ber of 3 ± 2 of the 250 \im sources were randomly associated with 
GNS galaxies in this test (where the uncertainty is the standard de- 
viation). This can be treated as an upper limit, as it assumes no 
correlation between objects detected in the sub-mm and near-IR - 
and so the real fraction of spurious matches is likely to be lower. 

Table [T] lists the properties (redshift, stellar mass, rest-frame 
colour) and flux densities of the individual detected sources. The 
median redshift of the detected objects is z = 2.02, and the me- 
dian stellar mass of the detections is log(M*/MQ) = 10.8. We note 
that in comparison to the bulk of the GNS sample (Section |2.1| ), 
these objects typically have lower photometric redshift probabili- 
ties, with median P = 61. 

Fig.[3]shows the location of the detected objects in the (U — B) 
colour-stellar mass plane. Clearly, relatively more massive galax- 
ies with red rest-frame (U — B) colours are detected, as shown 
in Fig. [4] We find that roughly 13 per cent of the sample with 
log(M*/M0) > 11 and (U -#) rest > 0.85 (the fiducial colour cri- 
terion adopted for dividing quiescent and star forming galaxies in 
Kriek et al. 2009) are detected at 250 \im. Given their far-infrared 
flux densities, these objects are clearly not quiescent, and we ex- 
pect them to have high dust masses and high star formation rates, 
with their red colours being as a result of dust extinction. How- 
ever, it is possible that the dominant origin of the IR emission is 
hot dust associated with AGN, rather than star formation, although 
this is not likely: e.g. Symeonidis et al. (2010) found that all of 
their 70 \im selected galaxy sample were primarily powered by star 
formation, Although X-ray AGN were removed from the sample 
at the outset (Section |2.1| ), we checked for additional AGN using 
colours in the Spitzer IRAC bands ( Stern et al. 20 05}, using dat a 
from the GOODS Spitzer Legacy program ^Dickinson et al.|2 003). 
Fig.[5]shows the [3.6]-[4.5], [5.8]-[8.0] colour-colour plot of the 
250 \im detected GNS galaxies. We find that six objects fall within 
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Figure 3. Distribution of 1.5 < z < 3 GNS galaxies with log(M*/MQ) > 
9.5 in the (U — B) rest colour-stellar mass plane (small red dots). The large 
black diamonds indicate the objects detected at 250 fim in HerMES. The 
typical uncertainty in the GNS stellar mass estimates is ~ 0.2 dex, while the 
typical uncertainty in (U — B) rest is 0.15 mag (see Conselice et al. 2011). 




Figure 4. Fraction of GNS galaxies with log(M*/M0) > 9.5 and 1.5 < z < 
3 detected with SNR > 3 at 250 as functions of rest-frame (U — B) rest 
colour (left) and stellar mass (right). Clearly, massive galaxies with redder 
colours are preferentially detected. For comparison, the rest-frame colour 
separation between quiescent and actively star forming galaxies adopted by 
|Krieketal.|j2009t is at (U-B) Test = 0.85 (dashed line). 



the region typically occupied by AGN. We do not remove these 
objects from the sample, as some studies have shown that AGN 



et al 



mainly contribute to the IR flux at wavelengths < 20 ^tm ( Netzer 
|2007[ |Mullaney et al.||2011| see also |Hatziminaoglou et al, 



2010 



we will instead note these objects in the following analysis 
(see also S ection [?. 3 . 3 | below) . 

We note that it is possible that the presence of either an AGN 
or starburst may lead to the stellar masses of some of the detected 
sources being overestimated. Other studies, which explicitly cor- 
rect for the effect of power law emission from AGN, find that ne- 
glecting such corrections can lead to differences of 10-25 per cent 
in stellar mass estimates of SMGs (e.g. Hainline et al. 2011). We 
show in Section |3.3.3| below that more sophisticated SED mod- 
elling, using rather different assumptions to those used in deriv- 
ing the GNS stellar masses, verifies that the 250 \im detected GNS 
galaxies are genuinely massive systems, (see also the discussion 
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Figure 5. IRAC colour-colour plot of GNS galaxies detected in HerMES. 
Overplotted are non-evolving tracks of various spectral templates as they 
are redshifted from z = to z = 2 (see legend; the crosses indicate the z = 
end of each track), taken from the library of Polletta et al. ( 2007). The 
colours of most of the objects are not consistent with those expected of 
Type I QSOs (shown by the shaded area marked 'AGN' in the legend), and 
are more similar to those expected of star forming galaxies at this redshift. 



concerning stellar mass estimates of AGN hosting GNS galaxies in 
|Bluck ctal.|20TT) . 



3.2 SED fitting 

To estimate Lir and SFR for the SPIRE detected GNS galaxies, we 
fit their far-IR SEDs using a modified blackbody (e.g. Hildebrand 
[T983l|Blain et al.|2003) of the form, 

S v =AvP£(v,7d ust ), (1) 

where Z?(v, T^ust) is the Planck function, A is the amplitude, and (3 
is the emissivity index (fixed to (3 = 1 .5). In addition, the Wien tail 
is replaced with a power law of the form S v oc v _oc , with a = — 2 
(Blain et al. 2 003]). W e also fit the SEDs using the templates of 
|Chary & E lbaz ( 2001 CE01 hereafter), as a consistency check on 
our results. 

We fit the SEDs using % 2 minimisation, allowing the dust tem- 
perature to vary in the range 10 — 70 K. We ignore the 24 \im flux 
densities when fitting the SEDs using models of the form of equa- 
tion [T] since at z > 1.5 we do not expect the modified blackbody 
model to be a reasonable description of the SED at this wavelength 
in the observed frame. However, we do include the 24 \im fluxes 
when fitting to the CE01 templates, as these include the contri- 
bution from poly aromatic hydrocarbon (PAH) features. Note that 
we include SED points with SNR< 3 in the fitting - given the re- 
quirement of a 24 \im detection and prior position, so long as the 
uncertainties on these points are accurately estimated, then the ad- 
ditional information they provide should help to better constrain 
the SED than either neglecting these points, or replacing them with 
3a upper limits. We comment on the effect of this on our results in 
SectionIP 

We derive the total (8 — 1000 |^m) IR luminosity (Lir) from 
the amplitude of the best fitting model, and convert this to a SFR, 
assuming that the Kennicutt (1998) law holds at this redshift, 

SFR IR (M yr- 1 ) = (4.5 x 10~ 44 ) x L lR (erg s" 1 ), (2) 



defined with respect to a Salpeter (1955) IMF. We therefore apply 
a correction of —0.23 dex to SFRs estimated using equation [2] to 
account for the Chabrier ( 2003) IMF assumed in this work (see e.g. 
|Kriek et al.|2009) . 

We also estimate dust masses during the SED fitting, using 
the method of Dunne et al. ( 2011 see also Dunne et al. 2000] and 
references therein), i.e., 



K 2 50^(v,r dlls 



(3) 



where 5250 is the A ux density at 250 \im in the observed frame, 
K is the ^-correction to rest-frame 250 fj,m, Dl is the luminosity 
distance, and K250 is the dust mass absorp tion coefficient, taken to 
be 0.89 m 2 kg -1 as in Dunne et al. ( 2011]). There are many caveats 



for the dust mass estimates obtained in this way, such as: the un- 
certainty in the value of K250; the fact that equation [3] can under- 
estimate the true dust mass due to the presence of warm dust in 
galaxies being neglected in the modified blackbody model (equa- 
tion^; and the large ^-correction to the redshift range of our study. 
Although the absolute values of M<j us t are highly uncertain, we use 
the relative values obtained by this method to give an indication of 
the relation of Md us t with M*, assuming that the dust properties are 
similar in galaxies of different stellar mass in our redshift range of 
interest (see Section [43] >. 

We estimate errors on the parameters derived from the SED 
fits using Monte-Carlo simulations. For each observed SED we 
generate 1000 random realisations, assuming that the errors on the 
fluxes are Gaussian. For objects with only photometric redshifts, 
we simultaneously randomise the redshift of the fitted model SED 
according to the scatter of a z = 0.06 measured by Griitzbauch et al.| 
P011| >. We adopt the 68.3 percentile range from the distribution of 
parameter values obtained from the random realisations as the cor- 
responding ±1g uncertainty. 



3.3 Results 

3.3.1 Star formation 

The SED fitting shows that the GNS galaxies individually detected 
in HerMES are ULIRGs, spanning the range 1 1 .9 < log Lir (Lq ) < 
12.9, with mean log Lir = 12.36 ± 0.05 L , where the quoted 
uncertainty is the standard error on the mean. We estimate to- 
tal SFRs for these galaxies under the assumption that this corre- 
sponds to the sum of the SFR derived from the far-IR SEDs and 
the UV-based unobscured SFR measurements from lBauer et al.l 
(2011). We find that the mean total SFR for these galaxies is 
SFRir + uv = 280 ±40 Mq yr -1 . Removing the six galaxies with 
IRAC colours consistent with AGN has no significant effect: with 
these objects excluded, we find SFRtr + tjv = 260 ±50 Mq yr -1 . 
This is a factor of > 2 larger than the mean UV- slope extinction 
corrected SFR estimates from Bauer et al. (2011) for these same 
galaxies, i.e. SFRuv,corr = 120 ±30 Mq yr -1 . We obtain results 
within < la of these values for all of these properties if we take into 
account the fraction of potential spurious matches (Section [TT] ) in 
a Monte-Carlo fashion. 

We checked the sensitivity of these estimates to the adopted 
sub-mm selection criteria. We find consistent results for the smaller 
sample of 8 galaxies detected with SNR > 5 at 250 |^m (mean 
log Lir = 12.39±0.09 L , mean SFRi R+uv = 290±60 M yr" 1 ), 
and for the sample of 14 galaxies detected at SNR > 3 at 350 \im 
(mean log L IR = 12.34 ± 0.07 L , mean SFRi R+uv = 260 ± 
40 Mq yr -1 ). We also checked the effect of including SED points 
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Table 1. Properties of 1 .5 < z < 3.0 GNS galaxies detected at 250 am with S/N > 3. Flux densities (S\) are in mJy, and only wavelengths in common between 
both GOODS-N and GOOD-S are shown. The error bars on photometric redshifts (we do not show error bars on objects with spectroscopic redshifts, marked 
with f ) and stellar mass estimates are statistical only, and the typical uncertainty in (U - #)rest is 0.15 mag (see Conselice et al. 201 1 for details). 
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logM* 
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1.10 
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0.63 
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20.8 ±4.3 
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8.4 ±2.6 


5.9±3.1 
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2.31 ±0.20 


10.50±0.13 


0.70 


0.234 ±0.005 


2.1 ±0.4 


4.3±0.6 


8.7±2.6 


8.4±3.2 


16.6±4.2 


6220 


1.72±0.16 


11.01±0.19 


1.21 
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7.0±1.5 


18.4±2.6 


15.2±3.2 


7.0±4.3 


7475 


1.61f 


11.23 ±0.06 


1.15 


0.191 ±0.004 


1.8±0.4 


4.6±0.8 


12.9±2.6 


7.1 =±=3.3 




7970 


2.54±0.21 


10.99±0.13 


1.05 


0.264 ±0.004 


1.7±0.5 


4.8=r= 1.1 


11.8±2.6 


10.3±3.2 
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* = IRAC colours of this object indicates AGN may be present (see Fig. 5) 

f = spectroscopic redshift (taken from the compilations by Barger et al. 2008 and Wuyts et a l.|2008) 



with SNR < 3 in the fits (see Section [3^2| ) - replacing them with 3a 
upper limits, we obtain mean log Ltr = 12.40 ± 0.05 L©, with cor- 
responding mean 

SFR IR+U v = 300 ±40 M yr _1 , for the whole 
sample of 22 galaxies. 

Dividing the sample by rest frame colour, we see no evi- 
dence for different IR properties for galaxies detected at 250 |^m 
with red or blue colours, although of course the sample is very 
small. We find mean logL TR = 12.33 ±0.09 (SFR TR+ uv = 270 ± 
60 Mq yr _1 ) for the 1 1 galaxies with (U — B) rQSt > 0.85, and mean 
logL IR = 12.34±0.06 (SFR IR+U v = 26O±4OM yr" 1 ) for the 11 
galaxies with (U — 5) rest < 0.85. 

We conclude that SFRi R+ uv is significantly higher than 
SFRuv,corr for our sample. | Wuyts et al.| ( |2011a| ) also found that 
SFRuv,corr is underestimated compared to SFRi R+ uv for galaxies 
with similar total star formation rates and redshifts to our sample. 
However, several other recent studies find the reverse situation. For 
example, Murphy et al. (201 1) observed a sample of 0.66 <z< 2.6 
24 am selected sources with additional 70 am photometry, and 
found that their measurements of SFRuv,corr are a factor of > 2 
her than SFRi R+ uv- They concluded that the dust corrections 



applied to their sample (from the Meurer et al. 



1999^ relation) were 



overestimated for many objects. Nordon et al. (2010 ) found similar 
results from a study using PACS observations of massive galaxies 
at 1.5 < z < 2.5 in GOODS-N, finding SFRuv,corr is overestimated 
by a factor of about 2 for galaxies with SFRtjv > 40 Mq yr _1 , 
assuming a Calzetti UV attenuation law (note however that | Wuyts] 
|et al.||2Qj"Ta| showed that this result may in part be driven by the 
relatively bright K s < 22 limit adopted in |Nordon et al.|2010| ). rBuat| 
|et al.| p010| > reached similar conclusions from a study of 250 am 
selected z < 1 galaxies from HerMES with UV photometry from 
the Galaxy Evolution Explorer (GALEX) satellite. At lower redshift 
(z < 0.35), Wijesinghe et al. (2011 ) found only a weak correlation 



with large scatter between the UV slope ((3) and Li R /Ltjv> which 
would also lead to overestimated SFRuv,corr- However, for UV se- 
lected samples (e.g. Lyman break galaxies; LBGs) which are not 
ULIRGs, dust corrections from the local Meu rer et al.| (l999) rela- 
tion appear to be valid at z ~ 2 (e.g. Overzier et al. 201 1 Reddy 
|et al.|2012| . Reasonable agreement between SFRuv,corr and SFR 
derived from stacked radio and 24 am observations is also seen up 
to z ~ 3 for LBGs ( |Magdis et al.|2010) . 

We expect large IR-derived SFRs for the galaxies we detect at 
250 \xm given their redshift and the 3a flux limit, which is ~ 9 mJy 
at 250 fim in the GOODS-N field. This leads to a large Malmquist 
bias (with some flux-boosting due to the low SNR) in comparison 
to the UV-derived SFRs, which reach to ~ 1 Mq yr _1 ( [Bauer et ah| 
|2011| ). Fig. [6] shows the SFRj R limit as a function of redshift for 
a modified blackbody model SED (equation [TJ with T^ ust = 35 K, 
normalised to a 250 \im flux density of 9 mJy. Highlighted in this 
plot are the SFRj R and SFRuv,corr values for the SPIRE-detected 
galaxies; and clearly in most cases SFRuv,corr is much lower than 
the fiducial SFRj R corresponding to the 250 \im flux limit. This 
makes the comparison between these two SFR measures for our 
sample difficult to interpret. There is one clear exception, where 
SFRuv,corr is roughly a factor of 3 larger than SFRj R - this is ID 
5918, which, from inspection of the ACS imaging, seems to be a 
multiple component merger system, with regions of significant un- 
obscured star formation (see Fig. [7}. It may be that only one com- 
ponent of this system is the source of the FIR emission, but it is not 
possible to determine which using the current data. 

Fig. [8] shows the comparison of SFRi R+ uv and M* for the 
SPIRE-detected GNS galaxies with the wider GNS sample, where 
for the latter SFRtjv,coit is used as the estimate of the total SFR. 
We see that almost all of the SPIRE detected galaxies scatter above 
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Figure 6. Comparison of SFRir estimated for GNS galaxies detected at 
250 [im (black diamonds) with the approximate 3a flux limit as a func- 
tion of redshift (blue line, estimated assuming a modified blackbody SED 
with Tdust = 35 K), and the extinction corrected UV estimates (SFRuv,corr) 
for these same galaxies (cyan squares). The SFRuv,corr values of the entire 
GNS sample are plotted for comparison (small red dots). The dashed lines 



indicate corresponding 


SFR estimates for a given galaxy. 




ID = 5918 
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Figure 7. ACS (V, i, z) image (10" x 10") of the mulitple component sys- 
tem ID 5918 (left), the only galaxy in the sample with significantly larger 
SFRuv,corr than SFRjr of the GNS galaxies detected at 250 \xm (see Fig.[6j. 



the SFR-M* relation measured by Daddi et al. ('2007), which is as 
expected given the approximate SFRjr limit shown in Fig. [6] 

3.3.2 Dust properties 

For the 16 galaxies with flux measurements in all SPIRE bands, we 
find dust temperatures in the range 23-48 K, with mean 35 ± 6 K 
(where the quoted uncertainty is the standard deviation). Note how- 
ever that only 4 of these galaxies have SNR> 3 in all SPIRE bands, 
and so the individual temperature estimates are poorly constrained, 
with typical statistical uncertainty « 5 K. We find that replacing 
the SNR < 3 SED points in the fits with 3a upper limits (see Sec- 
tion |3.2| ) gives T^ ust values for individual galaxies in this subsample 
that agree within < la of the values obtained when the low SNR 
SED points are included. For a sample selected with SNR > 3 at 




9.5 10.0 10.5 11.0 11.5 12.0 



log(M,/M ) 

Figure 8. Relation between total SFR and M* for GNS galaxies. The large 
diamonds represent SPIRE detected galaxies; those highlighted in blue have 
IRAC colours consistent with AGN (see Fig.[5j. For these galaxies, the total 
SFR estimate that we use is SFRir + uv- The small red points represent the 
wider GNS sample; in this case, the total SFR estimate is SFRuv,corr- The 
dashed line is the SFRuv, C on— M* relation measured at z ~ 2 by Dad di et aT] 
(2007). Note that the error bars indicate statistical errors in SFR and M* 
only. 

350 |^m, we find mean r<j ust = 33 ± 7 K, while for a sample with 
SNR > 5 at 250 \im, we find mean r dust = 34 ± 7 K. The single 
GNS galaxy which is detected at SNR > 3 at 350 \xm but is not in 
our 250 \im selected sample (ID 283; see Section [3T| ) has a slightly 
lower dust temperature (T^ mt = 20 ± 5 K). 

The mean temperature we find is somewhat lower than the 
typical temperature of ULIRGs at z < 1 (Tdust ~ 42 K, e.g. Cle ments | 
|et al.|2010||Yang et al.|2007| ); although note that p is fitted for in 
the former work, whereas in the latter it is fixed at (3 = 1.5, as we 
assume here. This is not unexpected given the high redshift of the 
sample and the selection at SPIRE wavelengths ( Symeonidis et al.| 
|2011| . The dust temperatures we find are similar to those found for 
other samples at z > 1 (|Chapin et al.|2009[|Amblard et al.||2010| 
|Chapman et al.|2010||Hwang et al.|2010| ). Adopting (3 = 2.0 in the 
modified blackbody model (equation[TJ gives mean T^ ust about 4 K 
lower. 

We find a mean dust mass for these galaxies of M^ust ~ 
3 x 10 8 Mq, which is comparable to the characteristic mass in the 
dust mass function of Mj ust ^4x 10 8 Mq measured at z ~ 2.5 
by [Dunne, Eales & Edmunds ( 2003 note the value quoted here is 
taken from Table 3 of Dunne et al. 2011). However, the range in 
M^ust spans more than an order of magnitude, and the individual 
values are highly uncertain. The median M^ ust /M^ ratio for these 
galaxies is ~ 5 x 10 -3 and spans the range 4 x 10 -4 — 3 x 10 -2 . 
This is similar to the range found by Rowan-Robinson et al. ( 2008), 
with a sample reaching to z ~ 2 and using a different method to es- 
timate Md us t. Fixing the value of (3 = 2.0 in the modified blackbody 
model (equation]!]) would increase the mean dust mass that we find 
by « 60 per cent. 

3. 3. 3 Joint optical-IR SED fitting 

We tested the sensitivity of the results described above to the sim- 
ple modified blackbody model used in the SED fitting (Section [T2| 
by jointly fitting the full optical-IR SEDs (BVizH from HST, IRAC 
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Figure 9. Examples of optical-IR SEDs fitted with CI GALE. Note that dif- 
ferent underlying assumptions were used with CI GALE compared to the rest 
of this work; i.e. the Maraston ( 2005) stellar population models, Kroupa 
pOOT) IMF, and |Dale & Helou|p002) infrared templates were used. The 
CI GALE fit results suggest that the bulk of the IR emission is associated 
with star formation rather than AGN. Note that the median %j: ed of the sam- 
ple is 1.7, so the example fits we show here are representative, although we 
choose to show ID 4180 in particular because it is the object with the largest 
inferred AGN contribution to the IR luminosity. 



channels 1-4, MIPS 24 \im, plus the Herschel photometry) using 
CI GALE ( [Noll et al. ||20091 ), a code which fits the attenuated opti- 
cal light from stars and dust emission associated with star forma- 
tion and AGN simultaneously. The available models for use within 
CI GALE differ from those assumed for deriving the GNS stellar 
masses (see Section [2T| and the SFRs estimated in this work (see 
Section [T2| . We used the Maraston ( 2005} stellar population mod- 
els to fit the optical part of the spectrum, the Dale & Helou ( 2002) 
templates to fit the dust emission, and a Kroupa ( 2 001| ) IMF. Some 
example SED fits are shown in Fig. [9] 

We find that CI GALE gives stellar masses that span the range 
10.0 < log(M*/M ) < 11.5, with median log(M*/M ) = 10.9, 
confirming that these systems have high stellar masses, as mea- 
sured in the GNS using a different SED fitting code (Conselice 
|et al.|201l) . A two sample Kolmogorov-Smirnov (KS) test reveals 
that the stellar mass distributions are not significantly different 
(p = 0.33), although there is a scatter of 0.23 dex in the residuals 
between the two stellar mass estimates for each galaxy. 

The SFRs estimated by CI GALE are systematically lower than 
the results obtained using the modified blackbody model, presum- 
ably as a result of the different stellar population model, IMF, 
and dust emission spectral templates used, but the mean SFR 
(210 ±30 Mq yr -1 ) is similar to that found from the modified 
blackbody SED fits, despite this. The fraction of the IR luminosity 




Figure 10. Distribution of stellar masses with redshift for the GNS cata- 
logue in both GOODS fields. The blue dashed lines indicate the samples 
used in the stacking analysis presented in this paper. 



due to warm dust associated with AGN estimated by CI GALE spans 
the range 3 — 30 per cent, with median « 5 per cent. This suggests 
that star formation is the primary source of the IR emission in these 
objects, as found in other studies (e.g. |Netzer et al.[ 2007 ; Mullaney 
|etal.|2011) . 



4 STACKING 

As shown in Section [TT] only 2.5 per cent of the 1.5 < z < 3, 
logM*(MQ) > 9.5 galaxy sample is detected in the 250 |^m maps 
used in this work, and the detected galaxies are ULIRGs with large 
stellar masses (~ 10 11 Mq). We therefore performed a stacking 
analysis to extend our study to galaxies with lower stellar masses 
and fainter far-IR luminosities. An additional advantage of the 
stacking analysis is that the results are less biased than those ob- 
tained from a small number of sources detected at low SNR. The 
stacking was performed on maps from which sources were not sub- 
tracted. Note that in contrast to the analysis in Section]?] additional 
maps at longer wavelengths than SPIRE were used in the stacking 
analysis (see Section [2^2]). 



4.1 Sample definitions 

We divide the 1.5 < z < 3 GNS galaxy sample into four bins of 
stellar mass, reaching to the log(M*/MQ) > 9.5 limit to which the 
survey is complete (Grutzbauch et al. 2011 ; Mortlock e t al.|201l) . 
Fig. [To] shows the location of the mass-limited subsamples in the 
(M*, z) plane, compared to the full GNS catalogue covering both 
GOODS fields. Due to the low SNR of the resulting stacked detec- 
tions (see Section [43] >, we are not able to divide the sample into 
redshift bins, nor examine subsamples of passive versus actively 
star forming galaxies (although note that the latter is investigated 
using the GNS galaxy sample by |Bauer et al.|201 1\ using UV-based 
SFR measurements). Table[2]lists the properties of the mass-limited 
subsamples we stack. 
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Table 2. Properties of the mass limited galaxy samples for GOODS-North, GOODS-South, and the combined sample. N indicates the total number of galaxies 
that were stacked in each sample; N Zspec is the number of these objects with spectroscopic redshifts; (z) is the median redshift of the sample; Nx is the number 
of objects which are detected in X-rays (these are not included in the stacks and are not counted in N). 



Mass Sample 
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North 

^spec 


(z) 


N x 


N 


South 

^spec 


(z) 


N x 


N 
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9.5 <log(M*/M ) < 10.0 


275 


30 


2.4 





233 


20 


2.4 


5 


508 


50 


2.4 


5 


lO.O<log(M*/M ) < 10.5 


105 


23 


2.2 


5 


111 


17 


2.3 


6 


216 


40 


2.3 


11 


10.5 <log(M*/M ) < 11.0 
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Figure 11. Result of stacking on random positions for each stellar mass bin 
in the GOODS-N SPIRE maps. The dashed line in each subplot indicates 
the stacked mean flux recovered when stacking on the real object positions, 
as listed in Tabled 



i 



4.0 

2.0 

0.0 

-2.0 

-4.0 
4.0 

2.0 

0.0 

-2.0 

-4.0 
4.0 

2.0 

0.0 

-2.0 

-4.0 
4.0 

2.0 

0.0 

-2.0 

-4.0 



" log(MJ > 11.0 y 

.tt- -i 








" 10.5 < log(Mj < 11.0 








" 10.0 < log(Mj < 10.5 

♦ . — i — -* — I-—' 

~- 1 1 1 1 1 1 


" 9.5 < log(Mj < 10.0 

i i 





100 200 300 400 500 600 

Wavelength (^m) 

Figure 13. Difference between the stacked flux densities in GOODS-N and 
GOODS-S for each stellar mass bin. Within the large uncertainties there is 
no significant difference between the two fields, although the stacked flux 
densities are generally fainter in GOODS-S. 



4.2 Method 

The far-IR data used in this work has low angular resolution, par- 
ticularly in the SPIRE bands where the beam sizes are 18 /r , 25 /r , 
and 36" at 250, 350, and 500 \im, respectively, resulting in rela- 
tively large confusion noise. The source densities of GNS galaxies 
per beam are also large (median 9 sources per beam at 250 |^m), 
and if the effect of clustered confused sources is not accounted for, 
the resulting stacked fluxes will be biased. 

We use the global stacking and deblending algorithm of |Kur-| 
|czynski & Gaw iser (2010 KG2010 hereafter) to mitigate the ef- 
fect of this bias (for other approaches to this problem see |B ether- 1 
|min et al.||2012| |Bourne et al.||2012| ). We generalised the method 
to simultaneously stack and deblend all of the mass-limited sam- 
ples (see Table [2|, in addition to two 'non-target' samples of ob- 
jects. The first of these non-target galaxy samples is drawn from 
the 24 |^m catalogue of Magn elli et al.| ( |2009| which is also used 
to provide prior positions for source extraction in the Herschel 
maps used in this paper). This catalogue provides coverage out- 
side of the GNS footprint, and allows infrared bright galaxies be- 
yond the edges of the GNS fields to be deblended. Since 24 \im 
bright sources are correlated with sources detected at PACS and 



SPIRE wavelengths, these objects are the most likely to contami- 
nate stacked flux measurements of the mass-limited samples at far- 
IR wavelengths. The second non-target galaxy sample consists of 
all GNS galaxies which are not 24 \im sources and not included 
in the stellar mass selected samples (i.e. with z < 1.5 or z > 3, 
and/or log(M*/M0) < 9.5). X-ray detected objects that are not in- 
cluded in the stellar mass selected samples (Section \Ta) were also 
included in this second non-target sample. 

We estimate errors on the stacked fluxes by bootstrapping: we 
run the stacking and deblending algorithm 1000 times, assigning 
the flux at each object position uniformly at random (with replace- 
ment) from the observed fluxes in each sample. During this process, 
the positions of all sources in the samples are kept fixed, and so the 
attenuation factors used in deblending sources (cc^- in KG2010) re- 
main constant (i.e. it is only the flux values that are bootstrap resam- 
pled). We adopt the 68.3 percentile as the uncertainty in the stacked 
flux. We also estimated errors by jackknifing (i.e. from the distribu- 
tions of stacked fluxes obtained after removing a single source from 
each stacking sample in turn), finding slightly smaller error bars - 
the detection significances inferred using the jackknife error esti- 
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Figure 12. Recovery of Lir and Ld us t when applying the stacking algorithm and SED fitting on simple simulated maps. 



mates are 0.1— 0.2a higher than those obtained using the bootstrap 
error estimates. 

We test the robustness of the mean stacked flux measurements 
by randomising the object positions in each of the stacking samples 
(both target and non-target samples) and running the stacking algo- 
rithm, repeating this process 1000 times. For simplicity, we perform 
this test using the GOODS -N sample only. We show the results for 
each of the stellar mass samples in the SPIRE bands (since these 
are the most likely to suffer from the effects of confusion as they 
have the largest beams) in Fig. [TT] With the exception of the low- 
est stellar mass bin, we find that the probability of a chance spu- 
rious stacked detection is higher for the lower resolution channels. 
The detection probabilities inferred from this null test are consis- 
tent with those obtained from stacking on real object positions and 
assuming the bootstrap error estimate; the maximum difference is 
0.3a, with detection significances inferred from the random stack 
tests being higher. 



4.2.1 Simulations 

We perform simple simulations to check that we can recover SED 
parameters such as Lir and T^ ust without significant bias. We create 
simulated maps with the same pixel scales as the real GOODS-N 
maps and insert Gaussian sources with the appropriate FWHM for 
each channel at the positions of real objects in the GNS catalogue. 
The simulated sources are modelled using the modified blackbody 
SED (equation[T]). We note that this is somewhat idealised, as we do 
not include different SEDs from those used in the fitting procedure. 

For the stellar mass selected samples, anticipating the Lir 
measurements obtained for the real maps (shown in Section [43] ), 
we set each model SED to have log Lir(Lq) =11.0, 11.5, 11.7 
and 11.9 for galaxies in stellar mass bins log(M*/M0) 9.5-10.0, 
10.0-10.5, 10.5-11.0 and > 11, respectively. We draw L dust for 
each galaxy in each stellar mass subsample from a uniform dis- 
tribution, with a slightly different (^dmt "^dusf^ ran S e use d for each 
bin: (15-45 K), (20-50 K), (25-55 K), (30-60 K), in ascending or- 
der of stellar mass. This ensures that each bin has different mean 
^dusb f° r clarity in the right panel of Fig. [12] 

Models for galaxies in the non-target sample of 24 |^m bright 
sources have log Lir(Lq) = 11, which is the median value we find 
for these sources when estimating their Lir from their 24 \im flux 



densities alone (where we estimate Lir for each source as the me- 
dian value over the full range of CE01 templates). We do not in- 
clude the non-target galaxies that were not detected at 24 |^m in 
the simulated maps. Each model source is redshifted to its corre- 
sponding z in the GNS catalogue. We apply a Gaussian random 
scatter of (1 +z p )x0.06 in redshift to galaxies with only photomet- 
ric redshift estimates to simulate the effect of incorrect redshifts, 
where the amount of scatter is as found by Grutzba uch et al.| < [201 1| > 
from a comparison of a subset of GNS galaxies with spectroscopic 
redshifts (see Section [TT}. For sources in the 24 \im detected non- 
target sample without redshift information, we assign their model 
SED a redshift selected at random from the redshift distribution of 
GNS galaxies detected at 24 \im. 

Fig. [12] shows the results of running our stacking and SED 
fitting code (Section [T2| on the simulated maps. We find that we 
recover Lir to within ±30 per cent down to the lowest stellar mass 
bin. We see that there is a small positive bias in 7d ust , with the re- 
covered value being at most about 7 K lower than the mean input 
^dust- This bias is absent if we set r<j ust to a fixed value for all galax- 
ies in each bin, and is likely to be a consequence of the smearing 
of the stacked SED shape due to the different redshifts and dust 
temperatures of the model SEDs that go into each stack. 



4.3 Results 

Table[3]lists the mean stacked flux densities for each stellar mass se- 
lected subsample in each field. We find consistent results between 
the northern and southern fields given the large uncertainties, al- 
though the stacked fluxes in the south are typically fainter than in 
the north for most stellar mass samples (see Fig. [T3j. The stacked 
S/N values are low: in the north, we obtain « 2 — 3a detections 
across almost all SPIRE and PACS bands for only the two most 
massive stellar mass bins. However, the detection significance in- 
creases to « 4a in some channels for the second highest logM* bin 
when the combined sample is used. The SNR in the lowest mass 
bin is only « la across the PACS and SPIRE bands when using the 
combined sample. 

Despite the low SNR for each individual SED point, we pro- 
ceed to fit the SEDs, in order to derive rough estimates of Lir and 
SFRir for each stellar mass bin. We include the low SNR points in 
the fits, rather than excluding them, or treating them as upper lim- 
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Table 3. Stacked mean fluxes (in mJy) for 1.5 < z < 3 GOODS NICMOS Survey galaxies in stellar mass bins. Ellipses (...) indicate where the solution was 
negative and therefore unphysical. 



Sample: North 



Wavelength (pm) 


9.5 <log(M*) < 10.0 


10.0<log(M*) < 10.5 


10.5 <log(M*) < 11.0 


log(M + ) ^ 11.0 


24 


0(\Q _l o nnc 
yj.yjyjo nz l/.l/uj 


0.048 ±0.013 


0.056 ±0.011 


060 + 01 8 


100 


0.11 ±0.11 


0.26 ±0.20 


0.35±0.19 


R4 + 4^ 


160 


99 4-0 96 


0.71 ±0.51 


0.53 ±0.59 


i 87 4-0 86 


250 


0.32±0.39 


1.15 ±0.64 


2.66 ±0.89 


3.23 ±1.47 


350 


0.59 ±0.47 


1.20 ±0.77 


3.23±1.10 


3.18±1.39 


500 


0.68 ±0.40 


0.64 ±0.76 


1.85±0.94 


2.03 ±1.38 


1160 


04 + o 04 


0.03 ±0.07 


0.31 ±0.09 


0.41 ±0.17 






Sample: South 




Wavelength (pm) 


9.5 <log(M*) < 10.0 


10.0<log(M*) < 10.5 


10.5 <log(M*) < 11.0 


log(M + ) > 11.0 


24 




0.020 ±0.006 


0.065 ±0.013 


f)4Q _i_ n 09 s 
yj.yj'-ty zn l/.l/z,j 


/ u 




0.03 ±0.04 


0.39±0.16 


OzL 4- 10 


100 




0.06 ±0.08 


0.98 ±0.42 








0.25 ±0.34 


2.72 ±1.20 


Q9 -I- 70 


250 




0.66 ±0.49 


3.52±1.11 


1.67 ± 1.21 


350 




0.80 ±0.57 


3.72±1.21 


2.15±1.58 


500 


0.06 ±0.36 


0.28 ±0.58 


2.80±1.11 


0.62 ±1.36 


870 


0.04 ±0.07 




0.29 ±0.19 


0.34±0.32 






Sample: Combined 




Wavelength (pm) 


9.5 <log(M*) < 10.0 


10.0<log(M*) < 10.5 


10.5 <log(M*) < 11.0 


log(M*) > 11.0 


24 


0.005 ±0.003 


0.033 ±0.007 


0.060 ±0.009 


0.056 ±0.016 


70 




0.03 ±0.04 


0.37±0.16 


0.07 ±0.09 


100 


0.01 ±0.06 


0.16±0.10 


0.66 ±0.23 


0.59 ±0.26 


160 




0.48±0.31 


1.60±0.68 


1.55 ±0.60 


250 


0.04 ±0.25 


0.89 ±0.42 


3.01 ±0.74 


2.68 ±0.95 


350 


0.26 ±0.29 


0.95 ±0.48 


3.38±0.82 


2.88 ±1.05 


500 


0.40 ±0.26 


0.41 ±0.46 


2.17±0.70 


1.59 ±0.95 


870 


0.04 ±0.07 




0.27 ±0.17 


0.31 ±0.27 


1160 


0.03 ±0.04 


0.04 ±0.07 


0.33±0.10 


0.43 ±0.19 



its. Under the assumption that the estimated error bars are reason- 
able (note that here they are obtained in a consistent way across all 
wavelengths), this should not bias the fit. We fit the SEDs for each 
stack using nearly the same method that was used for the SPIRE 
detected galaxies (Section [3^2] ). We make one change to the fitting 
procedure in order to account for the wide redshift range covered 
by the galaxy sample: during the Monte-Carlo procedure used to 
estimate error bars on the fitted parameters (i.e. Ljr, 7d ust , the un- 
certainties of which feed through to SFRjr and Md ust ), we bootstrap 
sample the redshift applied to the model SEDs from the distribution 
of redshift values in each stellar mass bin. This approximately dou- 
bles the size of the uncertainties on Ljr and SFR in comparison to 
those obtained when the redshift is held fixed at the mean redshift 
of the galaxy sample. Fig. [14] presents the stacked SEDs and best 
fit results using the modified blackbody templates for the northern, 
southern and combined samples. 



4. 3. 1 Star formation 

We obtain estimates of SFRjr for each sample with typically a fac- 
tor of two uncertainty, despite the low SNR measurements in each 
individual band. We find that the difference between the stacked 
flux densities measured for the GOODS-N and GOODS-S fields 
(Fig. [13]) leads to lower SFRs for most stellar mass bins in the 
GOODS-S sample. However, there is little tension between the 
SFRs measured in each field: the largest discrepancy is between 



the highest stellar mass bins, but even in this case the difference in 
the SFR estimates is significant only at the < 2a level. We find that 
the SFRjr estimates obtained using the modified blackbody model 
and the CE01 templates are consistent. 

We estimate mean total SFRir + uv for the stacked samples 
by adding to each sample the mean UV-based estimate of unob- 
scured SFR from |Bauer et al.| ( [2011| ) for the same galaxies in each 
stellar mass bin. Fig. |15| shows the resulting comparison with the 
mean UV-slope extinction corrected estimates (SFRuv,corr) from 
|Bauer et ah] p011| ) for the same galaxies. We see a rough agree- 
ment between the two measurements given the large uncertainties, 
although while in GOODS-N SFRir + uv is higher than SFRuv,corr, 
the opposite is true in GOODS-S. Much of this difference comes 
from a factor ~ 2 difference in SFRuv,corr between the two fields, 
with SFR uv ,corr being higher in GOODS-S than GOODS-N. For 
all stellar mass bins apart from log(M*/MQ) > 11, the difference 
in SFRuv,corr between the fields is significant at the « 3a level. 
The difference in SFRir + uv between the fields is less significant, 
at most 1.6a. Also, in GOODS-S, the highest SFR is seen for the 
2nd most massive log M* bin, in both SFRjr + uv and SFRuv,com 
although neither of these SFR estimates are significantly different 
from those measured for the most massive log M* bin. 

We checked for differences between the GOODS-N and 
GOODS-S samples that could lead to these effects. It is not likely 
that they arise from different redshift distributions: a two sample 
KS test gives p = 0.21, i.e. the distributions are not significantly dif- 



Far-IR properties ofz ~ 2 GNS galaxies 13 



North 



N = 25 

L IR = (7.3 ± 3.1) x 
SFR = 74 ± 32 yr 1 
T dust = 34.8 ± 7.0 K 
X 2 = 0.80 



10.0 < log(M„/M ) < 10.5 



N = 105 

L IR = (2.6 ± 1.3) > 
SFR = 26 ± 13 M 
T dust = 41.1 ±8.8 
X 2 = 0.29 



10.5 < log(M„/M ) 
1.5 < z < 3.0 
- N = 47 

L IR = (4.4 ± 2.4) x : 
SFR = 45 ± 25 M„ yr 
" a = 30.3 ± 5.3 K 
= 0.95 




; L IR = (1.2 ± 0.7) x 10 1 
SFR = 12 ± 6 M s yr 
r T dust = 34.6 ± 9.8 K 
= 0.99 




< log(M t /M ) < 10.0 

< z < 3.0 
= 275 



Wavelength (//m) 



Wavelength (/^m) 



L| R = (2.3 ± 1.5) x 
SFR = 23 ± 15 yr 1 
T dust = 30.8 ± 7.5 K 
X 2 = 1-14 



L| R = (6.5 ± 3.3) x 
SFR = 65 ± 33 M„ yr 1 
T dust = 38.9 ± 6.5 K 
X 2 =2.41 





10.0 < log(M„/M„) < 10.5 




1.5 < 2 < 3.0 


> 10° 


N = 111 


L| R = (1.1 ± 0.7) x 10 11 L 




SFR = 12 ± 7 M yr 1 




T dust = 31.1 ± 10"0 K 


Den: 


X 2 = 0.57 




lux 






Wavelength (//m) 



Wavelength (/itn) 



Figure 14. The stacked far-IR/sub-mm SEDs as a function of stellar mass 
in GOODS-N (top), GOODS-S (middle) and for both fields combined (bot- 
tom). Solid lines indicate the best-fitting modified blackbody model to each 
SED. 



ferent. Another possibility is environmental effects: the GOODS-S 
field contains a galaxy overdensity at z = 1.6 ( |Kurk et aT| [2009) 
which lies within our redshift range. This structure is thought to be 
a forming cluster of galaxies, and so the denser environment on av- 
erage relative to the GOODS-N field may lead to a higher fraction 
of quiescent galaxies in GOODS-S, and therefore lower average 



SFR. However, we find that excising the region within 2 Mpc pro- 
jected radius of this structure makes no significant difference to the 
derived SFRs. It seems likely that the difference between the results 
for each field can be ascribed to the small area covered by the GNS. 



4.3.2 SFR-M* relation 

Fig.[l6]shows the SFRir + uv - ^* relation obtained for the stacks in 
each field. Similarly to |Bauer et al.| ( |201 l) , we see a shallower rela- 
tion compared to the SFR-M* relation of Daddi et al. ( 2007 ), who 
measured SFRuv,corr ^ ^* ' 9 for star forming galaxies at z ~ 2. This 
is not surprising, because the GNS sample is selected differently, 
purely by stellar mass, and therefore includes quiescent in addition 
to star forming galaxies (see also |Bauer et al.|20TT| . Furthermore, 
the different methods used to measure SFR are also subject to dif- 
ferent selection effects. Using weighted least squares regression, 
we find the relation: 



log SFR uv +ir(M yr 



: (0.39 ±0.12) log (M*/M ) 

+ (-2.5 ±1.2) (4) 



for the combined GOODS-N and GOODS-S fields. The fits ob- 
tained for the individual fields are indicated in Fig.[l6]and are con- 
sistent within the errors. 

A straightforward comparison of our results with other works 
is not possible due to differences in sample selection, stellar mass 
completeness and the wide redshift range used here. |Karim et"ah] 
P011| ) performed a stacking analysis in 1.4 GHz data using a 
3.6 [jm- selected sample of galaxies in the Cosmic Evolution Sur- 
vey field (COSMOS; ||Scoville et al.|2007| ); however, this survey suf- 
fers from incompleteness for log(M*/M0) < 10.4 at z > 1.5. At- 
tempting a rough comparison of our measurements with this work, 
we find that our SFRir + uv estimates are systematically lower, 
for similar stellar mass and redshift ranges. However, the discrep- 
ancy is only significant at the 2-3a level for our most massive 
bin (log(M*/M0) > 11, where we find SFR a factor of ~ 2 less 
than |Karim et al.|2011| ), and there is reasonable agreement for the 
10.5 < log(M*/M ) < 11 bin. |Kurczynski et al.|(2010| > studied 
a sample of star forming BzK galaxies (sBzKs) in the Extended 
Chandra Deep Field South, comparing several methods of measur- 
ing SFR using essentially the same stacking algorithm we used in 
this work. Their sample was not stellar-mass-selected, but we find 
that our SFRjr + uv estimates for log(M*/M0) > 10.5 galaxies are 
in good agreement with their measurements (obtained using IR data 
from MIPS, BLAST, and LESS) at the same redshift as our study, 
after accounting for the Salpeter (1955) IMF assumed in Kurczyn- 
|ski et al.| ( [20T0] ). 



4.3.3 Ratio of obscured to unobscured star formation and 
relation to stellar mass 

We plot the ratio of obscured to unobscured star formation 
(SFRir/SFRtjv) as a function of stellar mass for the stacked sam- 
ples in Fig. [T7] Since the uncertainties are large, this is not well 
constrained from our data. For both GOODS fields combined, we 
find the relation 



log(SFR IR /SFR uv ) 



:(O.69±O.19)log(M*/M ) 

+ (-6.7 ±2.0) (5) 



using weighted least squares regression. As for the SFR-M* re- 
lation, the fits for the individual fields are consistent within the 



14 M. Hilton et al 



North 



♦ log(M„/M ) > 11.0 

• 10.5 < log(.:/./M .) < li.o 
■ 10.0 < log(.:/./M .) < 10.5 
A 9.5 < log(M„/M ) < 10.0 



20 40 60 80 100 120 140 

SFR uv , corr (M yr l ) 



South 



log(M„/M ) > 11.0 

10.5 < log(M,/M ) < 11.0 

10.0 < log(M„/M G ) < 10.5 



20 40 60 80 100 120 140 

SFR uv , corr (M yr - 1 ) 



Combined 



♦ log(M t /M ) > 11.0 

• 10.5 < log(.:/./M.) < 11.0 
■ 10.0 < log(M,/M s ) < 10.5 
A 9.5 < log(M,/M G ) < 10.0 



20 40 60 80 100 120 140 

SFR UVcorr (M yr x ) 



Figure 15. Comparison of mean SFR in each stellar mass bin derived from stacking (SFRir + uv; this work) with the mean SFR derived from the UV- slope 
extinction corrected rest-frame UV flux (SFRuv,corr)- The latter uses measurements described in Bauer et al. ( 2011). We calculate the mean SFRuv,corr using 
the same galaxies as in the stellar mass bins used in the IR stacking analysis, after first scaling the |Bauer et al.|(201 1) values to a |Chabrier|(2003) IMF. Results 
are shown for each GOODS field separately, as well as the combined sample. 




Figure 16. The relation between SFRir + uv and M* for galaxies stacked in bins of stellar mass (black diamonds). The blue line shows a weighted least squares 
fit to the relation. The dashed line shows the SFR-M* relation measured by Daddi et al. ( 2007) at z ~ 2 for comparison. The small red points show the 
SFRuv,corr measurements for individual GNS galaxies from |Bauer et al.ipoil) , highlighting the large scatter in this relation. The results are shown for each 
GOODS field separately, as well as the combined sample. 




log(M„/M ) log(M,/M ) log(M,/M ) 



Figure 17. The ratio of obscured to unobscured star formation (SFRir/SFRuv) as a function of M* for galaxies stacked in bins of stellar mass (where 
SFRuv is taken from the measurements of Bauer et al. 2011). The blue line shows a weighted least squares fit to the relation. The dashed line indicates 
SFRir/SFRuv = 1. The results are shown for each GOODS field separately, as well as the combined sample. 
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large uncertainties. The slope of this relation suggests that galax- 
ies with larger stellar masses on average have a larger fraction 
of obscured star formation compared to lower mass galaxies. A 
similar result is reported and discussed in Wuyt s et al.| p011b| ), 
who suggest that the mass-metallicity relation is responsible, with 
higher mass (metallicity) galaxies having larger dust column densi- 



ties and correspondingly larger SFRjr/SFRuv ratios (see also Pan- 



jnella et al.||2009| >. For galaxies with log(M*/MQ) > 11, we find 
the range spanned across the GOODS-N and GOODS-S fields is 
SFR IR /SFR UV - 6 - 20. For comparison, [Reddy et al.| ( |2012| ) find 
SFRir/SFRuv — 4.2 ± 0.6 for a sample of z ~ 2 Ljj V galaxies ob- 
served as part of the GOODS-Herschel project. 



4.3.4 Dust properties 

Although we derive estimates for T^ ust in each stellar mass bin from 
the SED fits (Section [3^ , they are not well constrained, with un- 
certainties ~10 K. All of the stellar mass samples in each field have 
Tdust m the 20-40 K range, consistent within errors across the stel- 
lar mass range, and consistent with the mean value found for the 
individually detected sources (Section [3.3.2| ). We note that simula- 
tions suggest that the T^ ust estimates from the stacked SEDs may 
be biased low, perhaps by roughly 7 K (Section [4.2.1| ). 

The estimates of Md us t we obtain are fairly low in com- 
parison to Mj ust , the characteristic mass in the dust mass func- 
tion, as m easured by|Dunne et al.|(|2011|) fo r < z < 0.5 and at 
z ~ 2.5 by |Dunne, Eales & Edmunds ( |2Q03| >. The largest value of 
M^t that we measure (~ 1.3 x 10 8 Mq), corresponding to the 
log(M*/M0) > 11 bin, is a factor of > 3 lower than Mj ust mea- 
sured bylDunne et al. (12003 1 at similar z, and also lower than M^ ust 
measured at 0.4 < z < 0.5 ( Dunne et al.|201l| ). This may be as a 
result of the purely stellar mass based sample selection used here, 
which contains both passive and actively star forming galaxies; nat- 
urally, the samples used in |Dunne et al . (2003 2011) consist of 
galaxies selected in the sub-mm, and are therefore dominated by 
dusty star forming galaxies. 

The relation we see between Md us t and M* is very poorly 
constrained (log Md us t oc log M^ 45±0 - 37 ), owing to the large un- 
certainties in the dust masses, but suggests a mildly decreasing 
Mdust /M* ratio with increasing stellar mass, with Md ust /M* falling 
from ~5x 10 -3 to ~ 7 x 10 -4 over the stellar mass range 9.5 < 



log(M*/M )<ll. 



5 CONCLUSIONS 

We have investigated the far-IR properties of a stellar mass selected 
sample of 1.5 < z < 3 galaxies drawn from the GOODS NICMOS 
Survey - the deepest //-band HST survey of its type prior to the 
installation of the WFC3 instrument - using deep Herschel 70- 
500 |^m photometry from the HerMES and PEP key projects. We 
found: 

(i) Only 22 galaxies from the sample are detected at SNR > 3 
at 250 \xm. They are ULIRGs (median log Ljr(Lq) = 12.4), have 
high stellar masses (median log(M*/MQ) = 10.8), and are located 
at z ~ 2. 

(ii) From fitting the SEDs of the SPIRE detected galaxies, we 
find they have mean SFRjr + uv a factor of > 2 higher than the UV- 
slope extinction corrected estimates of Bauer et al. (2011). How- 
ever, we note that the IR-based SFR estimate suffers from a signifi- 
cant Malmquist bias, making the interpretation difficult. The mean 



dust temperature of the 16 objects with flux estimates in all Her- 
MES and PEP bands (r d u S t = 35 ± 6 K) is slightly lower than found 
for ULIRGs at z < 1. 

(iii) Using a stacking algorithm which attempts to deblend 
sources, we find marginal detections (2 — 4a) at SPIRE wave- 
lengths when stacking the galaxy sample in bins of stellar mass, 
even for the highest stellar mass bins (log(M*/MQ) > 10.5). 

(iv) Despite the low S/N of the stacked flux measurements in 
each band, we obtain estimates of SFRjr for the stacked samples 
with factor ~ 2 uncertainties. We find that SFRjr + uv measured for 
the stacked samples is in reasonable agreement with measurements 
ofSFRuv corr for the same galaxy sample bylBauer et al.|(|2011|>. 



(v) We find a relatively shallow slope for the SFR- M* relation 
m2 -4±01 ) compared to previous studies (e.g 



Daddi et al. 



(SFR 

|2007| ), which is likely due to selection effects, as our purely stellar 
mass selected sample contains a mixture of passive and actively 
star forming galaxies. 

(vi) We find evidence for an increase in the ratio of ob- 
scured to unobscured star formation with increasing stellar mass 
(SFR IR /SFRuv oc M°- 7±0 - 2 ). This is most likely a consequence 
of the mass-metallicity relation, with higher mass and metallicity 
galaxies being more obscured. 

Since the far-IR and sub-mm data used in this paper are amongst 
the deepest that will be obtained by Herschel, it is clear that to 
make further progress in characterising the far-IR properties of low 
stellar mass (log(M*/M0) < 10) galaxies at z ~ 2 using Herschel, 
a much larger galaxy sample is needed, as will be provided by the 
Cosmic Assembly Near-infrared Deep Extragalactic Legacy Sur- 
vey (C ANDELS ; |Grogin et al.|2011[|Koekemoer et al.|20TT) . 



ACKNOWLEDGMENTS 

We thank the referee for many helpful comments which have im- 
proved this paper. We thank Amanda Bauer for providing the 
UV-based SFR measurements of GNS galaxies and useful discus- 
sions. MH and CJC acknowledge financial support from the Lever- 
hulme Trust and STFC. Support for the GNS was also provided by 
NASA/STScI grant HST-GOl 1082. 

SPIRE has been developed by a consortium of institutes led 
by Cardiff Univ. (UK) and including: Univ. Lethbridge (Canada); 
NAOC (China); CEA, LAM (France); IFSI, Univ. Padua (Italy); 
IAC (Spain); Stockholm Observatory (Sweden); Imperial College 
London, RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); and 
Caltech, JPL, NHSC, Univ. Colorado (USA). This development 
has been supported by national funding agencies: CSA (Canada); 
NAOC (China); CEA, CNES, CNRS (France); ASI (Italy); MCINN 
(Spain); SNSB (Sweden); STFC, UKSA (UK); and NASA (USA). 

PACS has been developed by a consortium of institutes led by 
MPE (Germany) and including: UVIE (Austria); KU Leuven, CSL, 
IMEC (Belgium); CEA, LAM (France); MPIA (Germany); INAF- 
IFSI/OAA/OAP/OAT, LENS, SISSA (Italy); and IAC (Spain). This 
development has been supported by the funding agencies BMVIT 
(Austria), ESA-PRODEX (Belgium), CEA/CNES (France), DLR 
(Germany), ASI/INAF (Italy) and CICYT/MCYT (Spain). 



REFERENCES 

Alexander D. M., Bauer F. E., Brandt W. N., Schneider D. P., 
Hornschemeier A. E., Vignali C, Barger A. J., Broos P. S., 



16 M. Hilton et al 



Cowie L. L., Garmire G. P., Townsley L. K., Bautz M. W., Char- 

tas G., Sargent W. L. W, 2003, AJ, 126, 539 
Amblard A., et al., 2010, A&A, 518, L9 
Bai L., et al., 2007, ApJ, 664, 181 

Baldry I. K., Balogh M. L., Bower R. G., Glazebrook K., Nichol 

R. C, Bamford S. P., Budavari T., 2006, MNRAS, 373, 469 
Barger A. J., Cowie L. L., Wang W.-H., 2008, ApJ, 689, 687 
Bauer A. E., Conselice C. J., Perez-Gonzalez P. G., Grutzbauch 
R., Bluck A. F. L., Buitrago R, Mortlock A., 2011, MNRAS, 
417, 289 

Bethermin M., Dole H., Lagache G., Le Borgne D., Penin A., 

2011, A&A, 529, A4 
Bethermin M., et al., 2012, A&A, 542, A58 
Blain A. W., Barnard V. E., Chapman S. C, 2003, MNRAS, 338, 

733 

Bluck A. F. L., Conselice C. J., Almaini O., Laird E. S., Nandra 

K., Grutzbauch R., 2011, MNRAS, 410, 1174 
Bolzonella M., Miralles J.-M., Pello R., 2000, A&A, 363, 476 
Bourne N., et al., 2012, MNRAS, 421, 3027 
Bouwens R. J., Illingworth G. D., Franx M., Chary R., Meurer 

G. R., Conselice C. J., Ford H., Giavalisco M., van Dokkum P., 

2009, ApJ, 705, 936 
Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Buat V., et al., 2010, MNRAS, 409, LI 

Calzetti D., Armus L., Bohlin R. C, Kinney A. L., Koornneef J., 

Storchi-Bergmann T., 2000, ApJ, 533, 682 
Caputi K. I., et al., 2007, ApJ, 660, 97 
Cava A., et al., 2010, MNRAS, 409, L19 
Chabrier G., 2003, PASP, 115, 763 

Chapin E. L., Pope A., Scott D., et al., 2009, MNRAS, 398, 1793 

Chapman S. C, Ivison R. J., Roseboom I. G., et al., 2010, MN- 
RAS, 409, L13 

Chary R., Elbaz D., 2001, ApJ, 556, 562 

Clements D. L., Dunne L., Eales S., 2010, MNRAS, 403, 274 

Conselice C. J., et al., 2011, MNRAS, 413, 80 

Daddi E., et al., 2007, ApJ, 670, 156 

Dale D. A., Helou G., 2002, ApJ, 576, 159 

Dickinson M., Giavalisco M., GOODS Team 2003, in R. Bender 
& A. Renzini ed., The Mass of Galaxies at Low and High Red- 
shift The Great Observatories Origins Deep Survey, p. 324 

Dunne L., Eales S., Edmunds M., Ivison R., Alexander P., 
Clements D. L., 2000, MNRAS, 315, 115 

Dunne L., Eales S. A., Edmunds M. G., 2003, MNRAS, 341, 589 

Dunne L., et al., 2011, MNRAS, 417, 1510 

Elbaz D., et al., 2010, A&A, 518, L29 

Elbaz D., et al., 2011, A&A, 533, A119 

Giavalisco M., et al., 2004, ApJ, 600, L93 

Griffin M. J., et al., 2010, A&A, 518, L3 

Grogin N. A., et al., 2011, ApJS, 197, 35 

Grutzbauch R., Chuter R. W., Conselice C. J., Bauer A. E., Bluck 

A. F. L., Buitrago R, Mortlock A., 201 1, MNRAS, 412, 2361 
Grutzbauch R., Conselice C. J., Bauer A. E., Bluck A. P. L., 

Chuter R. W., Buitrago P., Mortlock A., Weinzirl T., Jogee S., 

2011, MNRAS, 418, 938 
Hainline L. J., Blain A. W., Smail I., Alexander D. M., Armus L., 

Chapman S. C, Ivison R. J., 2011, ApJ, 740, 96 
Hatziminaoglou E., et al., 2010, A&A, 518, L33 
Hildebrand R. H., 1983, QJRAS, 24, 267 
Hwang H. S., et al., 2010, MNRAS, 409, 75 
Karim A., Schinnerer E., Martmez-Sansigre A., Sargent M. T., 

van der Wei A., Rix H.-W, Ilbert O., Smolcic V., Carilli C, Pan- 



nella M., Koekemoer A. M., Bell E. P., Salvato M., 2011, ApJ, 
730, 61 

Kennicutt Jr. R. C, 1998, ARA&A, 36, 189 

Koekemoer A. M., et al., 2011, ApJS, 197, 36 

Kriek M., van Dokkum P. G., Pranx M., Illingworth G. D., Magee 

D. K., 2009, ApJ, 705, L71 
Kroupa P., 2001, MNRAS, 322, 231 

Kurczynski P., et al., 2010, ApJ submitted (arXiv: 1010.0290) 

Kurczynski P., Gawiser E., 2010, AJ, 139, 1592 

Kurk J., et al., 2009, A&A, 504, 331 

Le Ploc'h E., et al., 2005, ApJ, 632, 169 

Luo B., et al., 2008, ApJS, 179, 19 

Lutz D., et al., 2011, A&A, 532, A90 

Magdis G. E., Elbaz D., Daddi E., Morrison G. E., Dickinson M., 
Rigopoulou D., Gobat R., Hwang H. S., 2010, ApJ, 714, 1740 

Magdis G. E., et al., 2011, A&A, 534, A15 

Magdis G. E., Rigopoulou D., Huang J.-S., Pazio G. G., 2010, 
MNRAS, 401, 1521 

Magnelli B., Elbaz D., Chary R. R., Dickinson M., Le Borgne D., 
Frayer D. T., Willmer C. N. A., 2011, A&A, 528, A35 

Magnelli B., et al., 2009, A&A, 496, 57 

Maraston C, 2005, MNRAS, 362, 799 

Meurer G. R., Heckman T. M., Calzetti D., 1999, ApJ, 521, 64 
Mortlock A., Conselice C. J., Bluck A. F. L., Bauer A. E., 

Grutzbauch R., Buitrago P., Ownsworth J., 201 1, MNRAS, 413, 

2845 

Mullaney J. R., Alexander D. M., Goulding A. D., Hickox R. C, 
2011, MNRAS, 414, 1082 

Murphy E. J., Chary R.-R., Alexander D. M., Dickinson M., Mag- 
nelli B., Morrison G., Pope A., Teplitz H. I., 2009, ApJ, 698, 
1380 

Murphy E. J., Chary R.-R., Dickinson M., Pope A., Prayer D. T., 

Lin L., 2011, ApJ, 732, 126 
Netzer H., Lutz D., Schweitzer M., Contursi A., Sturm E., Tacconi 

L. J., Veilleux S., Kim D.-C, Rupke D., Baker A. J., Dasyra K., 

Mazzarella J., Lord S., 2007, ApJ, 666, 806 
Nguyen H. T., et al., 2010, A&A, 518, L5 

Noll S., Burgarella D., Giovannoli E., Buat V., Marcillac D., 

Munoz-Mateos J. C, 2009, A&A, 507, 1793 
Nordon R., et al., 2010, A&A, 518, L24 
Nordon R., et al., 2012, ApJ, 745, 182 
Oliver S., et al., 2010, MNRAS, 405, 2279 
Oliver S. J., et al., 2012, MNRAS, in press 
Oliver S. J., Wang L., Smith A. J., et al., 2010, A&A, 518, L21 
Overzier R. A., et al., 2011, ApJ, 726, L7 
Pannella M., et al., 2009, ApJ, 698, LI 16 
Papovich C, et al., 2006, ApJ, 640, 92 

Papovich C, Rudnick G., Le Ploc'h E., van Dokkum P. G., Rieke 
G. H., Taylor E. N., Armus L., Gawiser E., Huang J., Marcillac 
D., Pranx M., 2007, ApJ, 668, 45 

Peng Y.-J., et al., 2010, ApJ, 721, 193 

Penner K., Pope A., Chapin E. L., Greve T. R., Bertoldi P., Brod- 
win M., Chary R.-R., Conselice C. J., Coppin K., Giavalisco M., 
Hughes D. H., Ivison R. J., Perera T., Scott D., Scott K., Wilson 
G., 2011, MNRAS, 410, 2749 

Perez-Gonzalez P. G., Rieke G. H., Egami E., Alonso-Herrero A., 
Dole H., Papovich C, Blaylock M., Jones J., Rieke M., Rigby J., 
Barmby P., Pazio G. G., Huang J., Martin C, 2005, ApJ, 630, 82 

Pilbratt G. L., Riedinger J. R., Passvogel T., Crone G., Doyle D., 
Gageur U., Heras A. M., Jewell C, Metcalfe L., Ott S., Schmidt 
M.,2010, A&A, 518, LI 

Poglitsch A., et al., 2010, A&A, 518, L2 



Far-IR properties ofz ~ 2 GNS galaxies 

Polletta M., et al., 2007, ApJ, 663, 81 
Reddy N., et al., 2012, ApJ, 744, 154 
Rodighiero G., et al., 2011, ApJ, 739, L40 
Roseboom L, et al., 2010, MNRAS, 409, 48 
Roseboom I. G., et al., 2012, MNRAS, 419, 2758 
Rowan-Robinson M., et al., 2008, MNRAS, 386, 697 
Salpeter E. E., 1955, ApJ, 121, 161 
Scoville N., et al., 2007, ApJS, 172, 1 

Stern D., Eisenhardt R, Gorjian V., Kochanek C. S., Caldwell N., 

Eisenstein D., Brodwin M., Brown M. J. I., Cool R., Dey A., 

Green P., Jannuzi B. T., Murray S. S., Pahre M. A., Willner S. P., 

2005, ApJ, 631, 163 
Swinyard B. M., et al., 2010, A&A, 518, L4 
Symeonidis M., et al., 2010, MNRAS, 403, 1474 
Symeonidis M., Page M. J., Seymour N., 2011, MNRAS, 411, 

983 

Viero M. P., et al., 2012, MNRAS, 421, 2161 
WeiB A., et al., 2009, ApJ, 707, 1201 
Wijesinghe D. B., et al., 2011, MNRAS, 415, 1002 
Wuyts S., et al., 2011a, ApJ, 738, 106 
Wuyts S., et al., 2011b, ApJ, 742, 96 

Wuyts S., Labbe I., Schreiber N. M. E, Franx M., Rudnick G., 

Brammer G. B., van Dokkum P. G., 2008, ApJ, 682, 985 
Yan H., Dickinson M., Eisenhardt P. R. M., Ferguson H. C, Gro- 

gin N. A., Paolillo M., Chary R.-R., Casertano S., Stern D., 

Reach W. T., Moustakas L. A., Fall S. M., 2004, ApJ, 616, 63 
Yang M., Greve T. R., Dowell C. D., Borys C, 2007, ApJ, 660, 

1198 



